Strained graphene in the presence of external probes 
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We revise tlie tiglit binding approach to strained or curved graphene in the presence of external 
probes such as Photoemission or Scanning Tunneling Microscopy experiments. We show that extra 
terms arise in the continuum limit of the tight binding Hamiltonian which can not be accounted 
for by changes in the hopping parameters due to lattice deformations. These material independent 
extra coupHngs are of the same order of magnitude as the standard ones and have a geometric origin. 
Moreover, we show that no /3-independent pseudomagnetic fields exist in strained graphene. The 
physical relevance of the extra terms to actual experiments involving local probes is discussed. 
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One of the most interesting aspects of graphene is the 
tight relation between its morphological and electronic 
properties. Although this issue has been explored at 
length in the theoretical literature [TJ [2], and there is 
a fair amount of related experiments [SHS] , recent results 
[TDHH] have given an extra push to the subject that will 
be explored in this work. 

In the continuum limit of the standard tight binding 
(TB) approach, lattice deformations couple to the elec- 
tronic excitations in the form of fictitious gauge fields or 
scalar potentials [11[T3|. The parameter that links the TB 
approach with the continuum elasticity theory, /3, is re- 
lated to the electron-phonon coupling and appears in the 
definition of the strain-induced fictitious fields. (5 reflects 
the changes in the hopping parameter t of the TB model 
with the changes of the relative distances between atomic 
nearest neighbors due to the lattice deformations. The 
continuum TB model has been very successful and has 
given rise to strain engineering in graphene [TH - fTH] . The 
experimental observation of "pseudo Landau levels" asso- 
ciated to specific strain patterns [7j [101 has been modeled 
with the TB pseudomagnetic fields and constitutes a very 
beautiful experimental confirmation of the theory. 

In a recent work 1191 inspired by a geometric approach 
to curved graphene [20n24] . the continuum TB Hamil- 
tonian was supplemented with additional /3-dependent 
terms arising from a higher order derivative expansion, 
which can be interpreted as a position-dependent Fermi 
velocity and an extra vector field. The presence of /3 in 
the coefficients of the strain-dependent terms is a con- 
sequence of the basic assumption of the discrete nearest 
neighbor TB approach, where the only parameter in the 
model is the hopping integral tij among nearest neigh- 



bors. The issue of the possible relevance of the actual 
atomic positions to the TB Hamiltonian was raised in a 
recent publication '25' , with the conclusion that an extra 
/3-independent gauge field was needed to complete the 
model. 

In this work we will complete the analysis of the con- 
tinuum limit of the TB Hamiltonian that models the 
coupling of lattice deformations to electronic excitations 
in graphene, with special emphasis on the observational 
consequences for current experiments. We will derive ad- 
ditional /3-independent terms needed to supplement the 
TB Hamiltonian in order to make a quantitative com- 
parison with the experimental data when the continuum 
approximation is a good description. In particular we 
will show that, although the additional terms modify 
the coefficients of the position-dependent Fermi velocity 
and extra vector field obtained in [TU] , no /3-independent 
pseudomagnetic fields [53] exist in strained graphene. 

In what follows we will assume for simplicity that there 
are no short range interactions or disorder connecting 
the two Fermi points, so that the low energy description 
around each point remains valid. As is well known in the 
TB-elasticity approach [Ull^, elastic deformations of the 
lattice give rise, in the continuum limit, to pseudovector 
potentials that mimic the coupling of real magnetic fields 
to the electronic current. The standard TB Hamiltonian 
in the continuum limit is 

Hrs^~^vo|d^.^^aM+^A,)^. (1) 

where vq = 3 to/2 is the Fermi velocity for the perfect 
lattice, with t the hopping parameter for nearest neigh- 
bors and a the lattice constant; j=l,2 (summation over 
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a repeated index is understood over the article), and aj 
are the Pauh matrices. The pseudogauge potential Aj is 
related to the strain tensor by 
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^yy 



) , A2 = ^i-'2u,y), (2) 



where f3= |91ogi/91oga|. The strain tensor is defined as 
Uij = I {diUj + djUi + dihdjh), where Ui and h are in- 
and out-of-plane displacements respectively. Note that 
one usually assumes that crystal deformations are small 
and ([ij is valid only up to 0{ufj) corrections. We will 
follow this practice for the rest of the paper. 

The prefix "pseudo" reflects the fact that, unlike real 
magnetic fields, pseudomagnetic fields of elastic origin are 
invariant under time reversal and couple with opposite 
signs to the two Fermi points of the graphene system. 
In both cases, observable effects are not associated to 
the vector potential itself, but to the field B = V x A. 
These are the vector fields giving rise to oscillations in the 
density of states discussed in the experiments [101 E] and 
in the various works on strain engineering. 

As shown in [TH] , if one uses the TB approach to go one 
order higher in the derivative expansion, the Hamiltonian 
([T]) becomes 



H' 



TB 



i / (fxtp''[vij{x)aidj+vo(Jiri+ivocriAi]'il;, (3) 



where the field Ai is the one given in ([2|, 
tensorial and space dependent Fermi velocity. 



Vij is the 



given by 



, and Ti is a new vector field 



dj Uij 



(4) 



The key observation in this paper is that TB Hamil- 
tonians describing strained graphene , and ^ in 

particular, are commonly derived in a specific reference 
system, the "crystal frame". The reason is that the 
Bloch waves = J^x &~^^'^cLx used to diagonalize the 
TB hamiltonian are written using the atomic equilibrium 
positions {a;}, which are regularly spaced and indepen- 
dent of the crystal deformation. On the other hand, in 
the presence of strain the positions measured in the "lab 
frame" are the actual positions of the atoms fji. The two 
sets of coordinates are related by yi = Xi + Ui{x), where 
Ui is the in-plane horizontal displacement vector. Note 
that the vertical displacements h are identical in both 
systems. In the classical theory of elasticity, crystal (lab) 
frame coordinates are usually referred to as Lagrangian 
(Eulerian) coordinates PT". 

Thus, the TB hamiltonian ([s]) is actually the crystal 
frame hamiltonian Hc{x). In order to describe the in- 
teraction of electrons with external probes or fields, we 
must use the lab frame hamiltonian HLab{y), i-e-, the TB 



Hamiltonian has to be rewritten in lab frame coordinates. 
The TB hamiltonian is the sum of the Dirac hamilto- 
nian Hq plus the terms induced by the lattice deforma- 
tions. As these are already 0{uij), we have to compute 
change-of-frame corrections only for the M^j-independent 
piece {Hq)c of the crystal hamiltonian. The computation 
is simplified by using the symmetric convention for the 
derivatives of the fermion fields 



(iJo)c = -ivo / d xipl{x)a^ di ipdx) 



(5) 



where diip = l/2{tp^diip — (9i?/'^)-0 and the subscript 
in ipc indicates that this is the fermion field operator in 
the crystal frame. The derivatives transform according 
to 



d _ dyk d 
dxi dxi dyk 



((5jfc + d^Uk)dk {5ik + u,k + t^£ik)dk, 

(6) 



where Uik — {diUk + dkUi)/2 is the linear piece of the 
strain tensor and LVEik = {diUk — dkUi)/2. We also have 
to transform the integration measure 



d'x = 



det 



dxi 



dyi 
(1 - u.,i)d^y 



d^y = |det((5ifc - Uik ~ ujeik)\d^y 



(7) 



On the other hand, "i/'l^'c is the particle density opera- 
tor in the crystal frame. As the number of fermions in 
any region should be frame independent, we must im- 
pose ipl'ipcd'^x = il'^ipd'^y, where ip{y) is the lab frame 
field operator. This implies 



1pc{x) 



det 



dxk 
dyt 



-1/2 



(8) 



which exactly cancels the Jacobian in ([t]) . The net result 



IS 



Wo / d'^xtpl{x)(Tiy^tpc{x) 



-ivo I d^y 'ip\y)ai^ip{y) + {uki + u!eki){ip^ ak^ip) 



(9) 

where the derivatives in the last term act only on the 
fermion fields. Finally, the dependence on the antisym- 
metric piece LoSij may be eliminated by a local rotation 
of the peudospinors 

V'(y) ^ e-'""^V'(y) ^ V'(y) - '^^^si^iy)- (lO) 

Indeed, the identity ia^a^ — Ekio'i shows that this rota- 
tion cancels the term proportional to in (|9|. A contri- 
bution proportional to 9fca; vanishes as well due to the 
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anticommutation relation {(T3, ak} 
yields 



for fc = 1,2. This 



HLab — HtB + Hceom 

where Htb is given by ([3| and 



(11) 



Geom 



ivQ / cPx^^ 



UkiUkdi + -{diUki)(Jk 



(12) 



In the last line we have used integration by parts to revert 
to the asymmetric derivative convention. Note that, to 
first order in the strain, /3-dependent terms are the same 
in both frames. Eqs. (Ill and (12 1 are the main results 
in this paper. 

As /3 ~ 2, the new /3-independent terms in Hoeom are 
of the same order of magnitude as those in the standard 
TB hamiltonian In particular, the space-dependent 
Fermi velocity derived in the TB formalism in [19' will 
become 



(2u, 



(13) 



with the corresponding correction for the vector field 



1 



djUrj. (14) 



Note, however, that these are not the only possible 
"frame effects" or consequences of the change of coor- 



dinates. Eqs. (11) and (12) give the Hamiltonian of iso- 



lated strained graphene. In order to couple the system to 
external fields such as electromagnetic waves in the con- 
tinuum limit, one has to use the lab frame coordinates. 
Thus the substitution di — > di 



in the lab frame hamiltonian (11). This gives an addi 



ieAi^^ has to be made 



tional strain-dependent but /3-independent correction to 
the hamiltonian describing the interaction with an elec- 
tromagnetic wave in the Coulomb gauge 

5H,,u = -Woe / (fxukiA1''\i;^cfki}). (15) 



The hamiltonian H^ab can also be obtained by per- 
forming the TB calculation directly in the lab frame. 
This derivation is explicitly given in the Supplemental 
Material, where a seemingly new additional pseudogauge 
field is shown to have zero curl. As is the only "new" 
vector field in strained or curved graphene, in what fol- 
lows we will comment briefly on its physical significance 
and compare it with the well known pseudogauge field 
Ai in Eq. ([2]). First of all, note that, unlike Ai^ Ti is 
not a functionally independent field. The reason is that 



the hamiltonian ^ is hermitian only for F^ 



—dv ■ 



Thus, a position dependent Fermi velocity requires the 
existence of the new vector field F^. 

A look at^j3| suggests that F^ is some sort of purely 
imaginary [28] counterpart to A^. However, this is obvi- 
ously wrong, as gauge potentials have to be real (hermi- 
tian) . The true nature of F^ is made apparent if we use 
the identity iukcr^ — Skicri to rewrite the relevant term as 
-ivocriVi = Woo-iFi, with 



Fi = F2CT3 , F2 = — Fi(T3. 



(16) 



Note that F^ is matrix- valued and hermitian. This shows 
that the vector field F^ plays the role a hermitian con- 
nection for the SO (2) group of local pseudospin rota- 
tions (10) generated by (T3. As a consequence, a posi- 



tion dependent Fermi velocity will be accompanied by 
pseudospin rotation ("pseudospin precession"), i.e., by 
electronic transitions between the two sublattices. This 
fact is even more obvious in the covariant model in |19j . 
where F.^ appears as the spin connection associated to 
fermions propagating in a curved background. Thus F^ 
is not a gauge field and can not give rise to the charac- 
teristic Landau levels of real or pseudo-magnetic fields. 
The only pseudogauge field in strained graphene is the 
well known Ai given by ([2]). 




w=(x^-y^,2xy) 



q;0.5 
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FIG. 1. Pictorial view of the strain field discussed in the text 
and the changes it produces in the density of states. The 
dotted (green) line represents the contribution from /3 depen- 
dent terms alone, while the thick (blue) line represents the 
total correction including frame effects. The three plots cor- 
respond to X — —L, 0, L for the displacement discussed in the 

text with Umax =0.2. 



To see the practical implications of this work to ac- 
tual measurements we now work out some practical ex- 
amples. Consider first a density of states measurement. 
The frame effects discussed are rather trivial in this case 
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but enough to exemplify the issue. The effect of the co- 
ordinate change wiU affect STM measurements when the 
tip resolution is large in units of the lattice constant (no 
atomic resolution). The local density of states (LDOS) in 
the lab frame can be computed approximately in the local 
limit, for a sufficiently smooth Uij. To do this, p{E,Uij) 
is computed assuming Uij is constant, and then its depen- 
dence on the position is restored in the final expression 
p{E,x) = p{E,Uij(x)). The LDOS can be computed in 
momentum space 



p{E) = j dqxdqytr{E - H{q^, qy))' 



(17) 



with the Hamiltonian (11) by diagonalizing Vij, which 

triE-Ho)-\ (18) 



amounts to a change of integration variables 

dq^dq^ 



p{E) 



with Hq the unperturbed Hamiltonian and v± the veloc- 
ity eigenvalues. This yields 



(19) 



(the factor of 4 is due spin and valley degeneracy) which 
to first order in strain can be computed to give 



p{E,x) = poiE) [(1-f /3tru-tr u)] 



(20) 



A simple but interesting example is provided by in- 
plane strains that are quadratic in the position, such as 
those associated to the triangular bumps that led to the 
observation of pseudo-Landau levels in STM [lU] and 
that have been explicitly produced in artificial graphene 
in [TT] . Remember that the TB gauge field associated to a 



strain tensor Uij is A cx (u; 



-2uxy). Consider first tion in (15) 



a deformation vector given by m = {x^ —y^, 2xy)umax/^L 
shown in the upper part of Fig. [T] It is easy to see that 
the associated pseudomagnetic field will be zero in this 
case. The trace of the strain tensor is tr u = Umaxx/L, 
hence a line scan along the y direction will give a per- 
fect constant V shape {p{E,x) ^ \E\), while along the 
x direction there will be a dilatation effect such that 
p{E,x) ~ (1 -f UmaxiP — '^)x/L)\E\, as depicted in the 
lower part of Fig. [T] for different values of x. Due to 
the frame effects discussed in this work there is an ad- 
ditional, material independent change in the magnitude 
of the LDOS that adds on top of the (3 dependent con- 
tributions. This is important to consider if one wants to 
measure the space-dependent Fermi velocity from a local 
probe with resolution larger than the lattice constant. 

An interesting thing happens if we now consider the 
same deformation vector but exchange Ux and Uy, i.e., 
u (X {2xy, — y'^). In this case there will be no volume 
effect (tr u=0) and the strain will give rise to a con- 
stant pseudomagnetic field whose associated density of 



states will show similar Landau levels oscillations along 
any scanline. A 90 degree rotation of the strain defor- 
mation will give the same V shape with a Fermi veloc- 
ity increasing this time along x=const. Finally, for the 
strain u oc {x^ — y^, —2xy) both the trace and the pseudo- 
magnetic field will be zero and there will be no effect al- 
together. It can be shown that the geometric vector field 
coming from the frame change does not affect the DOS 
at the linear order in Uij considered in this work. 

On the other hand, these examples are a simple illus- 
tration of the fact that the Honeycomb lattice is very 
anisotropic and, of course, does not have full rotational 
symmetry; hence similar looking deformations give rise 
to very different effects in the STM images. The impor- 
tant point is that, in the case of general strain, the frame 
effects discussed in this work will be responsible for ad- 
ditional spatial modulation of the intensity of the LDOS 
while preserving its energy dependence. 

Frame effects will also be important when the absolute 
orientation of the lattice changes locally. An example 
of this effect can be observed in the polarization depen- 
dence of ARPES signal [29]. The usual ARPES pictures 
of Dirac cones see only one half of the cones, due to the 
form of the matrix element of the lattice electron at the 
K point with the free electron that comes out. This effect 
sees the absolute orientation of the lattice: if the lattice 
is rotated with respect to the polarization of hght, the 
part of the Dirac cone that is observed also rotates. As 
before, in order to describe the physics in the lab frame, 
vectors in the crystal frame have to be rotated to the 
lab frame. This is again a /^-independent contribution. 
Note, however, that the suppression of part of the ob- 
served Dirac cones in ARPES is due to the interference 
between photoelectrons emitted from the two sublattices 
and, as such, goes beyond the continuum limit interac- 
Effects of local lattice rotations in ARPES 
have been reported recently in [3^. The frame effects 
associated to lattice rotations could also be observed in 
optical experiments like those described in [5T] . 

As a summary, we have shown that the TB descrip- 
tion of deformed graphene in the continuum limit must 
be supplemented with geometric terms originating in the 
change of coordinates needed to describe interactions 
with external probes. These are of course always present 
in the experiments. The correct Hamiltonian to use 
when trying to fit experiments is Hiab = Htb + ElGeom- 
The new terms are material independent and different 
from the usual gauge fields arising from deformation 
induced changes in the hopping parameter. The ex- 
tra terms associated with these frame effects are of the 
same form as those already present in the complete TB 
Hamiltonian ([s]), but come with /3-independent coeffi- 
cients. Moreover, the only "new" vector field in strained 
graphene is the connection P^ (also present in the ge- 
ometric formalism [19)). which is compatible with the 
symmetry analysis in |32H35) and required by the her- 
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miticity of the hamiltonian whenever we have a position 
dependent Fermi velocity. As discussed above, it is not 
a gauge field and will not give rise to Landau levels in 
the density of states, although it may have other physical 
effects, such as pseudospin precession. No material inde- 
pendent gauge fields having physical consequences arise 
in strained graphene. The frame effects described in this 
work will be relevant to local experiments with resolution 
A 3> a, for which a continuum limit is appropriate. 
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SUPPLEMENTAL MATERIAL 



The tight binding derivation in the lab frame. — 

The extra terms to be added to the standard TB calcu- 
lation due to frame effects can also be obtained by redo- 
ing the TB calculation directly in the lab frame. That is, 
we consider the TB Hamiltonian H = — X]<ij> tija\bj + 
h.c, but now we map the labels to positions in the lab 
frame 

^ = -^E4V4+h.c., (21) 
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where y = x + u{x) with x = miai + 771202, and Sn 
are the three nearest neighbour vectors (we foUow ref. 
[I9] for their definition and other conventions). As our 
interest here is in the /3-independent terms generated by 
the change of frames, we have assumed that the hopping 
parameters tij take their equilibrium value t. The real 



meaning of the relabeling in (21 ) is that non-equilibrium 



atomic positions are used in the Fourier expansions 



ife-[£+fi(a;)]^t 
fe' 



(22) 



Note that, due to the fact that crystal momentum k is 
purely two-dimensional, only the in-plane components u 
of a three-dimensional displacement (u(x), h{x)) will ap- 



pear in the Fourier expansions in ( 22 ). As a consequence, 
only the linear piece Uij of the strain tensor can give 
rise to frame effects, while the out of plane contribution 
dihdjh does not play any role in this regard. The same 
conclusion was reached in the main text by noting that 
only u enters the coordinate transformation that relates 
crystal and lab frames. 

To see how this will change the effective theory at 
the K-point, it is instructive to analyze ay further be- 
fore computing the Hamiltonian. If we restrict the states 
to k = K + Sk with Sk < A, we get 



^iK ■ X iK ■ u{x) 



E 

Sk 



e e ^ 'afc 



(23) 



and, comparing with the corresponding expression — 
^i{K+5k)-x^^ in the crystal frame, we observe two new 
contributions. The first one is e*^""*^^), which we can 
factor outside the integral. This is a trivial phase fac- 
tor that can be reabsorbed into ay by a local gauge 
transformation and has no effect on the physics. As 
shown below, if we do not reabsorb this phase, it will 
show up in the effective theory as a new gauge field 
-A-i = di{KjUj) — {uij + ijJtij)Kj. But this gauge field 
has zero curl by construction and produces no pseudo- 
magnetic fields, even for position dependent strains. The 
second term e'*'^ '"^^") cannot be eliminated by a gauge 
transformation in this way, and will induce extra terms 
in the Hamiltonian which precisely correspond to those 
in Eq. ^ after the field rescaling ([s]) is performed. 
Back to the actual computation, plugging (22) into 



(21 1 gives 



x,n k,k' 

(24) 

It is convenient to use a symmetric parametrization for 
the momenta: fc — > fc + q/2,k' — > fc — q/2, which cor- 
responds to the symmetric derivative convention in ([5|. 
Expanding to linear order in u yields 



x,n k,q 



e e "fe+g/2°*:-9/2 ^ 



--q • {u{x) + u{x + 5n)) — ik ■ {u{x) — u{x + (5„)) 



-I- h.c. 
(25) 



which, in terms of the Fourier coefficients of the displace- 
ment u{x) — J2q e^'^'^u{q) can be rewritten as 



ff = -t^e^('=-«/^)-^"a^^/,6,_,/2X 
Siq) ~ '-u{q) ■ (q{l + e^^'"'^) + 2fc (1 - e^^-'-\ 



h.c. 
(26) 



Expanding around the K-point and performing the 
sums over n as usual yields the matrix element 

Hk,q = ^[S{q)aiki + iqiUj{q)at{Kj + kj)-iqiUi{q)ajkj], 

(27) 

where we have redefined Sk — >■ k. Replacing iqiUj{q) = 
Uij{q) + uj(q)eij we finally obtain 

Hk,q = vo[S{q)(Jiki + (uij +ujeij)aiKj 

+ {uij + ujei-i)(7^kj - Una J kj] . (28) 

The last two terms are precisely those obtained from the 
direct coordinate transformation (|6])-([7]) of the contin- 
uum Dirac equation before the field rescaling ([S]), which 
eliminates the term proportional to uu. As anticipated, 
there is also a seemingly new gauge field 



(29) 



which is a total derivative and has zero associated mag- 
netic field p6]. It can be eliminated by a local gauge 



transformation -0 — > e 



-iK-u{x) 



